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Abstract 

T— I ■ 

, Quark and gluon distributions in the light-cone wavefunction of a 

I high energy hadron or nucleus are calculated in the saturation regime. 

■ One loop calculations are performed explicitly using the equivalence 

. between the parton distribution in the light-cone wavefunction and the 

I production distribution of that parton in a current-nucleon (nucleus) 

0> ' scattering. We argue that, except for some overall numerical factors, 

(—] ' the Weizsacker- Williams wavefunction correctly gives the physics of 

Q-i! the gluon distribution in a light-cone wavefunction. 

'f: ; 1 Introduction 

> . 

^ . The idea of parton saturation in QCD is at the heart of the interest in 

^ ! small-x hadron and nuclear physics. This idea has its simplest and most 

■ " " ' intuitive statement in terms of the light-cone wavefunction of a high energy 

hadron or nucleus. Saturation of quark and antiquark densities |^ is the 
statement that the density of quarks per unit area and per unit of two- 
dimensional transverse momentum, that is per unit of true transverse phase 
space, is limited by a constant times the number of colors so long as the 
quark momentum is below some momentum Qs, the saturation momentum. 
Above Qs the quark distribution becomes perturbative. The result is stated 
precisely in (29). While the value of Qs can depend on the particular hadron 
whose wavefunction is being considered and on the longitudinal momentum 

^This research is sponsored in part by the Department of Energy, Grant 
DE-FG02-94ER-40819. 
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fraction of the quark, the statement of saturation, as given in (29), contains 
no knowledge of QCD dynamics or of the hadron in question. The analogous 
statement of gluon saturation[0, 0, |^, ^ is somewhat different. Gluon densi- 
ties can be much larger than quark densities. The gluon density naturally has 
a term as a factor in its ultimate limit and a (inQs/(? also appears lead- 
ing to the expression (64) for the gluon density per unit of phase space in the 
light-cone wavefunction. Eqs.(29) and (64) are based on one-loop quark and 
gluon calculations, respectively, and for a large nucleus with an additional 
use of BFKL dynamics in the gluon case. However, because (29) and (64) 
are so directly related, in general, to the scattering of a quark-antiquark pair 
and a gluon pair, respectively, on a nucleus or on a high momentum hadron 
we believe that these results are quite general except for the possibility of a 
pure number, not depending on the hadron in question, as a multiplicative 
factor on the right-hand sides of (29) and (64). 

In the case of the gluon density in a large nucleus in the saturation region 
there is a very nice model suggested by McLerran and Venugopalan^ where 
the valence quarks of the nucleons of the nucleus are treated as the sources 
of Weizsacker- Williams gluons[§, |], ^] which make up the small-x gluon dis- 
tribuiton of the nucleus. This model leads to gluon saturation as expressed 
in (6) and (8). Except for an overall constant factor we believe these are 
general results and thus that the Weizsacker- Williams model is a good pic- 
ture of the gluon saturation regime of the high energy hadron or nucleus. In 
the Weizsacker- Williams approximation the saturated gluon distribution is a 
pure gauge field and this again is believed to be a general result 0. The fact 
that the saturated gluons are pure gauge fields does not mean that they do 
not have a direct physical interpretation. Indeed, our procedure for calcu- 
lating the gluon distribution in the nucleus is to note that it is the same as 
the spectrum of produced gluons[0, a physical object. However, the fact the 
Weizsacker- Williams gluons are pure gauge gluons is what allows an exact 
calculation of their distribution. 

The dynamics that leads to parton saturation is BFKL evolution ^ be- 
cause the increasing number of gluons that appear in a hadron's wavefunction 
due to longitudinal momentum or x-evolution occupy a common region in 
transverse phase space and so can naturally lead to the large values of 
(see (8)) required for saturation. It is possible that the turnover in 
observed recently at HERA|ll0|, [TI| at small x as is lowered below 2GeV^ 
may already indicate a saturation of the quark distribution [^, |r^, The 
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fact that diffractive production, 7* + P — > x + P, indicates an energy depen- 
dence much stronger than that suggested by soft physics |Tl|] may also be an 
indication that saturation has been reached at HERA[p!5|, |T6[. Deep inelastic 
scattering on nuclear targets at HERA should be a very good place to look 
for saturation. If the turnover in really is due to saturation one can 

expect that the turnover occur at > 5GeV^ in deep inelastic scattering 
off large nuclei. Another place where saturation effects may be important 
is in the very early stages of relativistic heavy ion collisions. At RHIC one 
expects the gluon saturation momentum to be about one GeV with minijets 
in that regime contributing most of the freed energy. At LHC the satura- 
tion momentum should be 2-3 GeV bringing saturation dynamics into play 
as a major determinant in the very early stages, well before equilibration, of 
heavy ion collisions. Finally, when potentials as large as A ~ are reached 
one has entered a whole new regime of nonperturbative QCD where, for ex- 
ample, instanton effects can become important. Since the large potentials in 
the wavefunctions are pure gauge fields perturbation theory remains valid in 
describing the light-cone wavefunctions]^]. However, in the very early stages 
of the central region of a head-on heavy ion collision these gauge fields are 
freed, and over the time of the freeing of the gluons large field strengths, 
Ffj,u ~ 1/(7, appear whose dynamics should be genuinely nonperturbative. 

In Sec. 2, we indicate how one can determine the quark and gluon distri- 
butions in the light-cone wavefunction by looking at quark and gluon pro- 
duction. The essence of the argument is that in a particular light-cone gauge, 
described in detail in Ref. [4], final state interactions are absent allowing par- 
ton production to be an indicator of the wavefunction of the hadron. In Sec. 2, 
the Weizsacker- Williams result for the wavefunction is briefly reviewed. 

In Sec. 3, we calculate the quark momentum distribution in a large nu- 
cleus. This calculation is equivalent to the one-quark-loop fluctuations in the 
Weizsacker- Williams background field of the nucleus. However, our method 
of doing the calculation, following Ref. [2], emphasizes the relationship be- 
tween the quark distribution of the nucleus and the cross section for scatter- 
ing a quark- ant iquark pair on the nucleus. Saturation then corresponds to 
blackness in the scattering of the quark- ant iquark pair on the nucleus. 

In Sec. 4, we determine the gluon distribution at the one-loop level. Here 
the saturated distribution includes a factor of inl/x replacing ^lnQl/(? from 
the Weizsacker- Williams result. 

In Sec. 5, we interpret the one-loop results in terms of unitarity limits and 
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the quantum mechanical shadow term from a black disc. 

Finally, in Sec. 6, we argue that higher corrections, beyond the one-loop 
correction should simply replace, up to a constant factor, the inVjx factor 
found in the one- loop calculation by -inQ^jl? leading again to the semiclas- 
sical result. 

2 Determining quark and gluon densities in 
the light-cone wavefunction 

The key observation allowing one to measure and calculate quark and gluon 
distributions in the light-cone wavefunction is that these distributions are 
directly related to quark and gluon production in hard scattering reactions 
initiated by currents coupling to quarks and gluons. It is a familiar result in 
deep inelastic electron-proton scattering that the structure function Fi gives 
a measure of the quark distributions 

/ 

at least in a first-order QCD formalism. We are now searching a stronger 
result. In (1) the transverse momentum of the struck quark is integrated over 
the range < £^ < Q^. What we now wish to determine is the quark and 
gluon distributions in a proton or nucleus for a definite value of the parton's 
transverse momentum, l^. These unintegrated quark and gluon distributions 
are determined by the cross section for producing quarks and gluons at a 
definite £_l in a deep inelastic reaction. The idea is to choose the boundary 
conditions, the ie's, of the light-cone gauge so that thre are no final state 
interactions [^] thus guaranteeing that the struck parton appears in the final 
state with unchanged momentum. 

2.1 The quasi-classical approximation 

This analysis has already been carried out in some detail]^ in a quasi-classical 
calculation of gluon production off a large nucleus in a deep inelastic reaction 
initiated by the "current" j = —\F^yF^u which couples directly to gluons. 
The result found for the spectrum of produced gluons, is most easily 
expressed in terms of the quantity 
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N{x) = I dHe-^^---f-^ (1) 



which is given asH, ^ 



!• ]\t2 1 

In (2) the saturation momentum Qg is given by 

Q's = ^^^^^^ P ^G{x, Q') (3) 

where p is the nuclear density and xG is the gluon distribution for a nucleon 
with 1/x^ = Q^, the scale at which gluons are measured, b is the impact 
parameter of the current-nucleus interaction while R is the radius of the 
nucleus. Eq.2 is valid in a quasi-classical limit in which Q'^-^xG{x, Q^) has 
neither Q^-dependence nor inl/x factors at small x. Two limits of (2) are 
noteworthy, (i) At small values of 

Nix) — >x2smaiia;GA(x,x^) (4) 

reflecting independent scattering on the various nucleons of the nucleus, (ii) 
For very large R 

NU)^n.r,.^^ (5) 

reflecting the saturation of the number of gluons per unit area. In momentum 
space, but neglecting the logarithmic dependence of Q^, 



This interpretation of saturation is made sharper by noting that 

-Jd%< A;^(&X(6 + X) >= nN{x) (7) 
so that (6) can be written as0] 

-<<(6M-te + ^)>=^i (8) 



showing that for gluons of transverse size Ax^ the maximum value of Ax^Aj^ 
is of order 1/a. This is gluon saturation. The <> in (7) and (8) indicates av- 
erages in the light-cone wavefunction of the nucleus. We note that satuation 
is different than shadowing [^] since (4) indicates that there is no shadowing 
in the quasi-classical approximation. Finally, comparing (1) and (7) we note 
that the number density of gluons in the light-cone wavefunction as given 
by (7) is equivalent to the distribution of produced gluons in deep inelastic 
scattering. 

In Ref. [4], detailed arguments were given that, with a proper choice 
of light-cone denominators, final state interactions are absent in light-cone 
gauge thus allowing gluon production, at a given transverse momentum, to 
directly reflect the transverse and longitudinal momentum distributions of 
gluons in the light-cone wavefunction. In Ref. [4], these arguments were given 
in the context of a quasi-classical (Weizsacker- Williams) approximation, how- 
ever the result appears to be more general as illustrated in Appendix A of 
this paper. In the next two sections of this paper we calculate first the pro- 
duced quark distribution and then the produced gluon distribution in deep 
inelastic scatterings off a nucleus in the one-loop approximation. The results 
can then be identified with the quark and gluon distributions in a nucleus in 
the one-loop approximation. 

3 Quark distributions 

In this section, we calculate the produced quark distribution in deep inelastic 
scattering off a large nucleus. Though the calculation is done at the one- 
loop level we shall argue that the result is quite general and valid also for 
scattering off protons at very small values of x. The calculation we are about 
to perform is not so far different from what has previously been donep| for 
the deep inelastic cross section. The new element which is added here is the 
determination of the transverse momentum of the leading quark which then 
gives the quark distribution in the light-cone wavefunction. 

3.1 The lowest order 

In order to set normalizations we begin by calculating deep inelastic scat- 
tering off a single nucleon and in the one-loop approximation. The relevant 
graphs are shown in Fig.l. We choose a frame where 
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Figure 1: Virtual Compton Scattering off a Nucleon in the One-quark loop 
Approximation. 

a 

and 

and, in addition, wc always suppose that g_ >> Q'^/2q_ so that the scatter- 
ing, say in a covariant gauge calculation, takes the form of the virtual (trans- 
verse) photon breaking up into a quark-antiquark pair which then scatters 
off the proton. The vertical lines running through the graphs of Fig.l indi- 
cate that the imaging part of the forward Compton amplitude is taken. It is 
straightforward to write 

xq + xq = j ^'^^[^^ + (1 - ^f] 



+ Q^z{i -z) {i + kf + g2^(i -z))^ [i^f dk^ ■ ^ ^ 

where 

6j = (6i,6V) = (0,0,e^) (11) 

is the polarization of the virtual photon while xG is the gluon distribution 
of the target nucleon. It is not difficult to see that (10) leads to the usual 
expression for the quark sea. In the logarithmic approximation k^ « ^ the 
term | P in (10) becomes 
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which after angular averaging in k becomes 

Again, in the logarithmic aproximation ^ « and z «\, with £^ of the 
same size as Q^ziX — z\ so that using (13) one finds from (10) 

x(g(x, Q') + g(x, Q^)) = £ ^^G(x, (14) 

which is the correct leading logarithnic form of the DGLAP[|l^, |18|, equa- 
tion thus confirming our normalization in (10). 

From (10) we an get the differential distribution in transverse momentum 



as 



dxq{x,Q ) <yQ ^--^ /"^ j r 2 , /1 \2i 

^ / i _ i + k \ ^<fk dxG 

■ [p + Q^z{l-z) ~ {l+ky + Q^z{l-z) J ' 1^ 'dB^' 

It is convenient to go from i to the conjugate coordinate x by using 

e-£ f (fx 

+ Q^z{l-z 

which gives 

, , i L + k ,,2 r (Pxid?X2 



(15) 



72 

/ 4^^' " ~ ■ V^o(VQ2x2z(1-z)) (16) 



+ Q^zil-z) {£ + ky + Q^zil- z)'' J 167r2 - " 

(1 - e-'^-^0(l - e*--^)e ■ v.,^o(V'Q'a;fz(l - z))e ■ v..^o(v'Q2a;iz(l - z))). 

(17) 

Thus, 

^ = ^1 - ^~"') + (1 - ^ - (1 - e-^^-(^^-^^))] 
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z^+{^-zf]-^.,MjQHlz{l - z))-^,MJQ^xlz{l - z))e 



(18) 

In the logarithnic approximation the exponential factors involving k in (18) 
can be expanded through second order with the resulting integration over k 
limited by the corresponding x— factor. Thus, 

:i-»-**)^|? = ^^^G(.,x;) (19) 

where x\G{x,x\) is an abbreviation for x\G{x,Q'^ = l/xf), leading to 

^ ]^^^4 j d^xid^X2dz[x\xG{x,x\)+xlxG{x,xl)-{xi-X2fxG{x,{xi-X2) 

■ [z" + (1 - zf] v., i^o ■ \7.,K,e-''--^^^-^-\ (20) 

3.2 The quark distribution for a large nucleus 

Now consider (20) as the single scattering approximation for quark produc- 
tion on a large nucleus. Introducing the nuclear density, p, and the impact 
parameter corresponding to the nucleon in the nucleus, 6, we may write (20) 

as 

dxq Q^Nc f ,2, ,2 j2 r^^l^ HT? T^Cf , Xjv r—, -t:C, 



d^i 647r6 



2A 



e-^-fei--2) . [^2 ^ ^1 _ ^fYz^.. Ko{^Q^xlz{l -z))- y,,,K,{^Q^xlz{l - z)). 

(21) 

In arriving at (21) we have used[^] 

^ = j^^P^l^G{^^^l) (22) 

where 

2 j d%pVR^ - = A (23) 



with A the atomic number of the nucleus. 

Now it is straightforward to allow the quark-antiquark pair coming from 
the virtual photon to scatter on an arbitrary number of nucleons in the 
nucleus. One simply makes the replacement Q 

Introducing, the saturation momentum, Qs, by 

. f v/fi^l (25) 



one finds 

dxq Q'^Nc 



rf2£ 647r6 

(26) 

We suppose » Q^. Then the dominant contribution to (24) comes from 
the region z << 1. It is convenient to define a scaled variable y = Q^z in 
terms of which 



dxq Nc 



I d^X^d'x2il + e-(^i-^2)'«'/4 - e-^?«'2/4 _ g-xiQ^/4^^_,,.(,^_,^) 



d%dH 647r6 

■ dy Vxi Ko i^/xly) ■ V^^a^o {^/xh )• (27) 

It appears difficult to give a closed form for all the integrals in (27). 
However, it is rather simple to evaluate (27) either when i"^ » Q"^ or when 
£^ << Qg. In these cases 



dV^ 67r4 P 
and 



In arriving at (29) we have used 



for t « Qi. (29) 



J d'x,d^X2e-''-<^^'^^Uy ■ X7.M^) = 16vr' (30) 
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and 



I d^x.ffx^e-^^^-^^^^'^'dy v., K,{^y) ■ V.M^) = ISvr' (31) 

while the other two terms in (27) are neghgible when £^ << Q^. The integral 
(31) coming from the second term on the right-hand side of (27) corresponds 
to the multiple scattering of the observed quark by the medium in both the 
amplitude and in the complex conjugate amplitude. The integral (30) coming 
from the first term on the right-hand side of (27) is the shadow term for (31) 
corresponding to complete blackness of the scattering of the quark- ant iquark 
pair on the nucleus in the region £^ << Q^. 

Eq.(29) reflects satuation of the quark density in the nucleus for £^ << Q^. 
There is, up to a constant, one quark per unit phase space in the saturation 
limit. Unfortunately, we are unable to give a physical interpretation of the 
constant appearing on the right-hand side of (28). As will be discussed 
in more detail for the gluon case we believe that, except for the constant 
factor, (29) is completely general and that it does not depend on the one-loop 
approximation which we have used or on the fact that we have considered a 
large nucleus rather than a hadron with the quark at a very small value of 

X. 

4 Gluon distributions 

Now we turn to determining the gluon distribution in our large nucleus. This 
has earlier been done in the quasi-classical approximation. Here, we do the 
calculation at the one-loop level, a calculation from which we shall be able 
to extract a general result in Sec. 6. The calculation we are about to perform 
is related to that done some time ago0, however, here we focus on the 
transverse momentum distribution of the leading gluon, a quantity related 
to the light-cone quantized distribution of gluons in the target. 

4.1 The lowest order 

In order to set normalizations we begin with the lowest order calculation. 
We use the "current" 

jix) = --f;Xu (32) 
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which produces gluons off a nucleoli at lowest order from the graphs illus- 
trated in Fig. 2. We parametrize q and £ as in (9) while the polarizations of 
the produced gluons are written as 



^1 <^ 



Ai 



'\2 

i2 



{L+k) 



(33) 



(34) 



and for shortness of notation we shall often write e'^^ = and e^"^ = e^- 
In what follows we shall always work in alogarithmic approximation for lon- 
gitudinal momentum so that we may assume £_ << g_. The sum of the 
graphs shown in Fig.2 have already been evaluated in Ref.[2] with result, 
when g2 » (2^ ^ ^Yien {i + kf « Q^z 





l+k 





(a) (b) (c) 

Figure 2: Lowest order graphs for gluon pair production off a nucleon. 



giving 



+ k 



If the graphs of Fig.2 are interpreted as light-cone perturbation theory graphs, 
in contrast to Feynman graphs, and evaluated in A_ = gauge then graph 
c is zero while the sum of graphs a and b is 



ra I -pb 

\, ii„ ~r i ) 



ei-{i + k)e2 -{i + k) e^-ie2- 



AiA2^ A1A2 {i + ky 

A simple calculation shows that 



£2 



(36) 
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Ai A2 

while 



(Z+IF ^^^^ 



E I^AiAa + TaiAzI^ - 2 £2^£^J^y (2^) 



A1A2 



SO that (37) and (38) are identical when k^/i^ « 1 after an angular average 
over directions of k. Thus in what follows we consider only graphs (a) + 
(b), which terms have an interpretation as the scattering of a color neutral 
two-gluon system on the target. 

For scattering of the current j off a single nucleon we again have four 
graphs exactly as in Fig.l but with j replacing the electromagnetic current 
jfj, and with the lines £,q — £,■ ■ ■ now referring to gluons rather than quarks. 
Then analogous to (10) one can write 

AaNr ^ r dH dz 



TT ^ ^ J 



e^- {l + kk2- il + k) .2d^k dxG{x,k'^) 



Ai,A2 ' 

2; 



(39) 



" e {i + kY " F dk^ 

where, again, we work in an approximation where there are no loops in 
xG{xik'^) so that k'^-^xG{x, k"^) is a constant, both in inl/x and in k"^. The 
Q^— dependence on the left-hand side of (39) comes from a cutoff < 
which is understood. Using (38) in (39) one finds 



xGix, Q') = r ^ xG{x, f) (40) 

which is correct in the leading double logarithmic limit, thus checking our 
normalization in (39). 

Going back to (39) one can write the differential distribution, the unin- 
tegrated gluon distribution, as 

dxG _ aN, ^ fiVQ' dz^e^-ie^-i ei • (£ + k)e2 -{l + k) ^^d'^k dxG 

2~" 

(41) 



r^yVc ^ /■« /"^ dz ei 
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We note that ^ is identical to what we ailed ^ in Sec. 2, and in Ref.[4]. 
Writing 

. • - = - / e-*--L(,, . - (42) 

one finds 

dxG _aNc 1 dz d?Xi(Px2 ^{xiX-iY ^ 

. (1 + f>-iHx^-^^) _ ^-ik-x, _ f,ik-x,^^l± (43) 
Duplicating the steps that led from (18) to (26) leads to 

dxG ^ Nl - 1 r dz d^xid^X2 (xi ■ _ 

dWe 87r6 ; 

(l + e-^^i-^2?Ql/^ _ e-^?Q?/4 _ e-^i«'/4) (44) 

where now 

Ql = ^^/W^ (45) 

which is identical to (25) except for the absence of the Cp/N^ factor and 
^ = §j with g as defined in Sec.2. 

Again, it is easy to evaluate (44) either when P » or when P « Q^. 
Thus, 



dxG Ql f dz - 1 _ g2 



dV£ 47r4 e 



J i = \^^nl/x^torf»Ql (46) 



When << only the 1 and e =^2)2Q^ ^g^j^s contribute to (44) exactly 
as happened in going from (27) to (29). Using 



and (see Appendix B) 
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/ ^^^^i^(2^^i^ - l)e-«?fe-^^)V4 = ^ (48) 

we arrive at 

c^xG _ Nl-l rQVQ' dz - 1 , , , . 

where the hmits of the 2;— integration are given by assuming the transverse 
momentum of each of the gluons approaching the nuclear target is of order 
Qs and by requiring the gluonic system have a coherence length > R. M 
is the nucluon mass. As in the fermion case the integral (48) corresponds 
to absorption of the two-gluon state as it passes over the nucleus while the 
contribution (47) can be viewed as the quantum mechanical shadow of that 
absorption and, as usual, the shadow and the absorption terms are equal 
when the target is completely absorptive (black). The new element here is 
the longitudinal momentum integral, the In 1/x factor in (49). In Sec. 6, 
we shall focus on what happens when the calculation is done beyond the 
one-loop level and what happens to the £nl/x factor in that case. 



5 Interpreting the one-loop results 

Now, however, let's try to understand the significance of the rather simple 
results contained in (46) and (49) as well as in (28) and (29). Eqs.(28) 
and (46) of course are straightforward and represent a process which is hard 
enough so that only a single nucleon in each nucleus is effective. Thus, (28) 
and (29) give quark and gluon number densities which are, after integrating 
over the impact parameter 6, just A times production off an isolated nucleon. 

It is the region where << which is more interesting. Refer for a 
moment to (27) where ^ is given as a sum of four terms on the right-hand 
side of that equation. The coordinate Xi refers to the transverse position 
of the observed quark in the amplitude while X2 refers to the same quantity 
in the complex conjugate amplitude. The second term on the right-hand 
side of (27), the e~^-i~-2) ^^/^ term, corresponds to the S-matrix for the 
quark- antiquark pair coming from the virtual photon and interacting with 
nucleons in the nucleus, both elastic and inelastic interactions, as it passes 
over the nucleus. Only the interactions with the observed quark do not 
cancel between real and virtual (inelastic and elastic) reactions. The fact 
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that {xi — ^2)^ < l/<5s shows that typically the measured quark will have 
transverse momentum on the order of Q^. Thus the contribution when i'^ << 
Ql is determined by the probability that a quark which gets many random 
"kicks" be found with relatively small transverse momentum. It is natural 
that dxq oc (Pi the phase space be available to the quark. Thus, except for 
normalization this is purely a statistical problem for £^ << Q^. The fact 
that ^ is independent of both Q"^ and Ql comes from the fact that the 
quarks which dominate the process are those having transverse momentum 
on the order of Qs before the quark-antiquark pair passes over the nucleus. 
Quarks having transverse momentum much greater than Qs, before reaching 
the nucleus, are not freed while passing over the nucleus while quarks having 
transverse momentum much less than Qg, before reaching the nucleus, are 
few in number and can be neglected. Since the total number of "effective" 
quarks is proportional to and distributed according to phase space the 
functional form, a constant, of (29) follows with one-half of that constant 
given by the second term on the right-hand side of (27). 

Still in the region £^ << Ql, the first term on the right-hand side of 
(27) corresponds to no scattering whatsoever of the quark-antiquark pair by 
the nucleons of the nucleus. It can be viewed as the quantum mechanical 
shadow of the term described just above. When a quark-antiquark pair 
having relative transverse momentum 2£ impinges on the nucleus, and if 
Ip « Ql, this pair always interacts with the nucleus with the momentum 
of the quark and antiquark getting distorted far from L The destruction of 
this part of the wavefunction is accompanied by a "shadow" term where the 
quark again has momentum L This is the first term on the right-hand side 
of (27). 

For the gluon-loop (44), (46) and (49) are direct analogies to (27), (28) 
and (29) with the new element being the longitudinal phase space of the 
gluons. When the "current" j breaks up into a gluon-gluon pair there is a 
large phase space for one of the gluons (the observed gluon) to carry almost 
all the current's longitudinal momentum while the other gluon carries a small 
amount which, however, because of the vector nature of the gluon gives a 
logarithmic integral in the probability that the current break into a gluon- 
gluon pair. It might seem that this gives an arbitrarily large factor as inl/x 
becomes large, but this is not quite so as we shall now see in the next section. 

Finally, we note that when < Ql and ^ have an interpretation as 
quark and gluon densities in a light-cone wavefunction but they do not have 
the interpretation as quark and gluon distributions in terms of an operator 
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product expansion. This is most easily seen in the discussion of Sec. 2 where 
^ = has significant A-depcndcnce, but where there is no shadowing 
whatsover and no nontrivial A-dependence in terms coming from the operator 
product expansion. 

6 What happens to the £n 1/x factor? 

In this section we shall argue that when higher quantum corrections are 
included the inl/x factor in (49) gets modified so as to lead to an expression 
essentially identical to (6). To see what happens to the in 1/x factor in 
(49) we must go beyond the simple small-x dynamics we have used so far 
in our discussion. We continue to find it useful to imagine the scattering of 
the current j on a large nucleus where we choose an unusual frame in order 
to explain, heuristically, what is the essential dynamics. Thus, choose the 
kinematics representing the system before the collision occurs, and illustrated 
in Fig. 3, to be 



q-k 




Figure 3: Kinematics of the current-nucleus scattering just before the colli- 
sion. 
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(50) 
(51) 
(52) 



with p the momentum per nucleon of the nucleus. With s = 2q_p_^_ fixed 
choose g_ large enough so that, for a given z, k-/k+ is significantly larger 

than one but not too large. That is, take g_ = A^^o ■ — with A^o a fixed, and 
moderately large, number. The idea here is to put most of the longitudinal 
momentum into p leaving just enough in q so that the process may be viewed 
as a two-gluon system, q — k and k, colliding with a highly evolved wave- 
function of the nucleus. By restricting the longitudinal momentum of the 
left-moving two-gluon system coming initially from j we need not consider 
further evolution in that system, at least in a leading logarithmic approxi- 
mation in longitudinal momenta. All logarithms except the dz/z integration, 
which is our focus, are included in the wavefunction of the nucleus. Now if z 
is decreased one must correspondingly increase g_ and decrease in order 
to keep g_ = N^k/ z and s = q-p+ fixed. By decreasing p+ we limit the range 
of useful x-evolution in the nucleus and in so doing decrease the saturation 
momentum Qg. But we must guarantee that k^ < in order that a reaction 
occur with reasonable probability, and this determines the lower limit of z 
in the dz/z integration. Our task then is to determine the x-dependence 
of Q^. In the quasi-classical approximation Q^, given by (3) or (45), has no 
x-dependence. However, once we go beyond the quasi- classical approxima- 
tion we expect an x-dependence. Indeed, (3) gives an x-dependence through 
xG{x, Q^) and it is to the determination of the x-dependence of this quantity 
that we now turn. 

For dynamics we use the fixed coupling BFKL equation which incorpo- 
rates leading logarithmic x-cvolution, including the small-x approximation 
to DGLAP evolution. Let xG{b,x,Q^) be the gluon number density for the 
nucleus at momentum fraction x, at scale Q, and at impact parameter b. Our 
normalization is such that 

J d% xG{b, X, Q^) = xGa{x, Q'^) (53) 

the normal gluon distribution of the nucleus. It is xG{b,x,Q'^) that we ex- 
pect to replace 2^JW^^ pxG{x,Q^) in (3) and (45). Then, in the BFKL 
approximation 

xG{b, X, Q') = a [ No{b, A, Q^)e2^^W^+^^- «V0^ ^ (54) 
where we expect A^o to be slowly varying in A. In Eq.(54) 
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X(A) = V'(l) - - li'il - A) (55) 

where the A-integration goes along the imaginary axis, while Y — i nl/x. 
The saddle point of the A-integration is determined by 



X (Ao) = 2aiVe Y ^^^^ 



givmg 



xG{h, X, Q') = (g)^V^-(^°)^. (57) 

We expect (57) to be valid so long as >> Ql with Ql being determined 
by 

.Gtt,.,® = £(|) (68) 

where xG{b, x, Ql), approached from the perturbative regime Q^ > Ql agrees 
with (49) when i'^ — > Ql from the lower momentum side, i'^ < Ql- We allow 
c to have weak (logarithnmic) x and Ql dependences. Using (57) and (58) 
one finds 

n \ \i) r^-iinfi o r ^^^^ \ , 2aiVcx(Ao) , . 

(1 - \)ln QJQ, = ^n{-j====} + Y (59) 

Using (56) in (58) gives 

The right-hand side of (60) is slowly varying in Y and in Aq. Thus, for very 
large Y, leading to very large Ql, Aq is determined by 

(ax) 

a value of Aq which is not too far from A = 1/x so that in 1/x evolution 
dominates Q^— evolution indicating that our appoach to the problem should 
be reasonable. 



19 



Turning to (59) and using the fact that the first term on the right-hand 
side of that equation is slowly varying in Y one can determine that 



1 dQljY) ^ 2aN,x{Xo) 
Ql dY 7r(l - Ao) 



(62) 



giving the dependence of the saturation momentum on Y. 

Now we are in a position to answer the question of what happens to the 
£n 1/x factor in (49) when higher quantum corrections are included. Let 
Y = i n s/Q^ and let Y{i) be that rapidity such Ql{Y{e)) = Then the 

integral in (49) becomes 



where en 1/ xq — Y — Y[€) so that, using (63) one gets (49) to become 



which, apart from the factor, is identical to (6). We feel that the essential 

factors in (64), the ^^^^ and the In Q^je? factors, are general results in QCD 
for the light-cone wavefunction. The overall constant in (64) we do not trust. 

It is perhaps useful to consider carefully why we claim that (49) is gen- 
eral, except for the issue of the In 1/x factor which we have discussed in 
some detail in this section. To that end turn to (44). Except for the four 
terms in parentheses at the end of the right-hand side of (44) all the other 
factors reflect the current j breaking up into a gluon-gluon pair along with 
the Fourier transform going from transverse coordinate to transverse mo- 
mentum space. Thus all the dynamics of the target is in the last factor. Of 
the four terms constituting the last factor the first, the 1, corresponds to no 
interactions of the gluon pair with the target and this term is universal and 
independent of the nature of the target. The third and fourth term corre- 
spond to possible interactions only in the amplitude and complex conjugate 
amplitudes, respectively. These two terms are functions only of x\ and x\, 
respectively and so must generally integrate to zero when £^ << Q^. The 
second term corresponding to S{xi)S*{x2), the product of the (S— matrices 
in the amplitude and complex conjugate amplitude, need not in general take 
the form given in (44). This term in general must be 1 when Xi = X2 and 
it should be small when {xi — X2yQl » 1, but we have no argument as 




(63) 



dxG _ Nl-l 1 - Aq 
dWe ~ ^M^c 2x(Ao) 



(64) 
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to the exponential form of (44) being exact. It is this specific form which 
gives (47) and (48) the same value. Physically, the exponential form in (44) 
came from the many independent scatterers in the nucleus. While this is 
also natural in our more general circumstance we do not know how to prove 
it. In any case, the expectation that the term in question go to zero when 
{xi —Xj^YQI >> 1 immediately gives a contribution which is £— independent 
when P « Q"^ since the £— dependence only comes in through e~*-'(-i~-2). 
The real issue then is what constant replaces the tt^ on the right-hand side of 
(48). We have written (64) as if that constant remains vr^ while it may pos- 
sibly be some other pure number. Thus the form given in (64) we feel must 
be true, but there may be an additional constant multiplying the right-hand 
side of that equation. 



Appendix A 

In this appendix we illustrate how, with a proper choice of ie's in light- 
cone denominators, final state interactons can be suppressed. The example 
given here is similar to that given some time ago ||2l| although now we have 
a better physical interpretation of what is happening. 

Consider the graph shown in Fig. 4 where there is a final state interac- 
tion of the gluon (k) with the struck quark. In light-cone gauge only the 

term pq— is important in the gluon propagator, p has a large + com- 
ponent of the momentum, and we suppose q has only -|- and — compo- 
nents with 2g+g_ = —Q^. We may assume that all lines in the upper blob 
have + components of their momentum greater or equal to (£ — — g)+ = 
a;p+while|(£-g)_| << g_. Thus {l-k-q^ ~ -{L-Jlf and P ^ -k^ so that 
one need only consider {l — kY + ie and /c+ as denominators possibly trapping 
the contour and thus limiting /c_|_ to small values. It is only in case the 
contour is trapped at a value where |fc+| < that final state interac- 
tions are important. Now (A; — £)^ + ie~ —2l_{{k — t)^ + {^^=^ — it). Thus, 
if the light-cone denominator is taken to be [kj^ — ze]^^ there is no trapping 
of the /c+ contour while any other choice leads to trapping. This choice of ie 
corresponds to the gauge potential extending to large negative x_— values, 
but not to large positive x_— values^, thus naturally avoiding final state 
interactions. It is straightforward to include additional gluons connecting to 
the struck quark. What is not so clear is whether or not there is a consistent 
definition of light-cone denominators which renders higher loops finite and 
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at the same time eliminates the gauge field for large positive values of x_ 
This is an important technical problem yet to be resolved. 




Figure 4: Potential final-state interaction which is absent with appropriate 
choice of boundary conditions for the light-cone gauge propagator. 

Appendix B 

In this appendix we outline how the integral in (48) can be evaluated. 



.2™2 V ^2™2 ' 1 



1-^2 



{Bl) 



1-^2 



then 



' - \ / ^ ^^^^ ^ ^^^4 ^^^^^ ^ i^^-^^i^ 

{B2) 

m 



Using 



1 r2n , 

27r io 



One finds 



I = 71 



0:1X2 2 



h{-QlxiX2)e 



-(xl+xl)Ql/i 



Now (See formula 19 on page 197, Ref. 



(54) 



j(°°^/2(ig^XiX2)e-^-^«^/^ = ^e^?«^/S(2,Q^x?/4) 
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where 7(a,x) is the incomplete 7— function. We find 



I = n' r ^7(2, z) = n'. (55) 
Jo 
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